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Abstract 

The recently observed excess in gamma-ray signal near the Galactic center suggests that dark 
matter particles may annihilate into charged fermions that produce gamma-ray to be observed. In 
this paper, we consider a leptonic dark matter, which annihilates into the standard model leptons, 
g + gT and r + r~, by the interaction of the gauged lepton number U(1) L _ L ^ and fits the observed 
excess. Interestingly, the necessary annihilation cross section for the observed gamma-ray flux 
provides a good fit to the value for the relic abundance of dark matter. We identify the preferred 
parameter space of the model after taking the existing experimental constraints from the precision 
measurements including the muon (g — 2), tau decay, neutrino trident production, dark matter 
direct detection, LHC, and LEP experiments. 
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I. INTRODUCTION 


The dark matter (DM) problem is one of the pressing issues in particle physics and cos¬ 
mology. While the existence of DM has been firmly established through various observations 
of its gravitational effects on multiple scales, its microscopic nature still remains unknown 
PP. This situation stimulates a variety of DM searches including the direct detection of 
dark matter scattering off detector materials, the detection of indirect signals from the dark 
matter annihilation or decay, and the collider searches of missing energy signatures due to 
the produced dark matter particles. Of particular, we notice that the new cosmic-ray detec¬ 
tion experiments, such as PAMELA [2], AMS-02 [3J, and Fermi-LAT [3], based on satellites 
reach unprecedented sensitivity to the cosmic-ray signals, which leads better chance to get 
the indirect information of dark matter properties. 

An intriguing observation was made using the public data of the Fermi Large Area Tele¬ 
scope (Fermi-LAT) by Hooper et. al. and also other independent groups [5TIT6] : a gamma- 
ray excess at ~ O(GeV) coming from the Galactic center (GC) is found. In the analyses, 
it is claimed that the gamma-ray excess spectrum is in good agreement with the emission ex¬ 
pected from DM annihilation into standard model (SM) charged particles^] The GeV excess 
is well fitted by a DM particle with a mass of tudm ~ 30 — 40 GeV annihilating to bb with an 
annihilation cross section of {av) ~ 2 x 10 -26 cm 3 /s [13, [16| |^] Silk et. al. pointed out that 
contributions of the diffuse photon emissions from primary and secondary electrons produced 
in DM annihilation processes are significant, especially for leptonic final states {it) [Hj. It 
is also noticed that with the inverse Compton scattering (ICS) and bremsstrahlung contri¬ 
butions from electrons, annihilations of DM particles with a mass of tudm ~ 10 GeV into 
ll provide a good fit with an annihilation cross section of {av) ~ (1 — 2) x 10 _26 cm 3 /s [13] , 
The bb final state may be understood by Higgs portal type DM models and studied by sev¬ 
eral authors [2T1 [21] but a model for the leptonic explanation based on leptophilic DM is 
relatively less studied for the GeV excess. Here we explore a leptophilic model with the DM 
mass ttt-dm ~ 10 GeV. 

In the heavier mass domain, Mdm ^ 100 GeV, leptophilic DM models have attracted a 

1 In Ref. 17], the authors proposed a new mechanism naturally inducing a continuum bump signature in 
cosmic 7 -ray measurements even with a particle directly decay into two photons, introducing Energy Peak 

idea together with the postulate of a generic dark sector m- 

2 We note that recent observation of AMS-02 [19] has started to exclude the \X bb dominant DM 


explanation of relic abundance [20] . 
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lot of attention (see e.g. [25].) due to recent observation of excessive cosmic-ray positron 
fraction by the PAMELA, Fermi-LAT, and AMS-02 experiments, but lack of excess in the 
anti-proton fraction |2b| . In building the leptophilic DM model, it is attractive to gauge 
the differences in lepton numbers: U(l)i e __L M , U(1 )x _z, T , and U(1 )l t -l b - These symmetries 
are anomaly free without extending the SM particle contents [271429] |^] Leptophilic DM 
models with a 11(1)^-^. gauge symmetry for the positron excess were studied in Refs. [501 - 
[32] . In our analysis, we take U(l)^ _^ T symmetry for the GeV gamma-ray excess since 
models with U(l)^ e _r , and U(l)£ T _£ e are stringently limited by existing cosmic-ray positron 
measurements in low energy [53] . 

It should be also noticed that astrophysical uncertainty in gamma-rays from around the 
GC including modeling of background emission in the inner galaxy is still big. Moreover, 
millisecond pulsars Ml Ho] [HUM] and pions from the collision of cosmic-rays with gas [Sh 
SHO] can contribute to the GeV scale gamma-ray and have been proposed as alternative 
explanations of the GeV gamma-ray excess even though the spectral shape from millisecond 
pulsars looks too soft at the sub-GeV energy range to account for the observed GeV excess 
spectrum [35]. Also the morphological feature of the observed excess is extended to more 
than ~ 10° from the GC beyond the boundary of the central stellar cluster which could 
include numerous millisecond pulsars [15], and observed distributions of gas seem to provide 
a poor fit to the spatial distribution of the signal HSIE6U371. 

The contents of the paper is as follows. In Section II, we explain the leptophilic DM 
model in detail and present dominant annihilation channels. The model parameter space 
for the observed DM thermal relic abundance and the GeV gamma-ray excess is clarified. 
In Section III, we discuss the existing experimental constrains on the same parameter space, 
then conclude in Section IV. 
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The other anomaly free choice is U(1 )b-l, but it does not provide lepton specific interactions. 
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II. THE MODEL, RELIC ABUNDANCE, AND FERMI-LAT GEV EXCESS 


A. The model 

We extend the SM: 

• by extending the gauge symmetry with U(l)r _ Lt , 

• by introducing a new Dirac fermion 0, which is identified as dark matter. 

The charge assignment for the SM fermions and the new fermion regarding the L^ — L r 
symmetry is given in Table [TJ The muon-leptons and anti-tau leptons are (+1), tau-leptons 
and anti-muon leptons are (—1) and the new fermion has a charge Q We take a universal 
gauge coupling constant g' for Z' interactions. 


particle 

0 

= (Vfj.L,gL), g-R, 

L 3 = (v t L,Tl), TR, V t R 

others 

charge 

% 

+1 

-1 

0 


TABLE I. Charges under the L^ — L r gauge symmetry. 

For the (spontaneously broken) extended gauge symmetry, we associate a new vector 
boson Z' with an undetermined mass my. The model Lagrangian is written as follows: 

C D Csm - + \m\,Z' a Z' a + 

+ g'QyZ'Jri’ + g'z' a V QpW, (i) 

where Q'^j are U(1 charges of the DM and a SM fermion /, respectively given in 

Table [T| In our study, the DM mass m^ is taken to be 10 GeV to fit the GeV excess as 
suggested in Ref. o (see also Ref. [38]). 

The 0 particle is neutral under the SM gauge interactions but its presence is seen by 
Lp — L r interactions. The gauge interaction allows an early time thermal equilibrium with 
the SM particles and the standard freeze-out took place at T ~ m^/ 20 through the the 
dominant annihilation channels: 

00 (—> Z'^) —> , Vfpn 

00 -> Z'Z 1 , 
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( 2 ) 

(3) 









where t = fi, r. The corresponding Feynman diagrams are depicted in Fig. [lj The DM 
annihilation into a Z' pair is kinematically allowed only when > mz>- 



£ + , v,. 


£ ,v e 



FIG. 1. Dominant DM annihilation channels: (a) s-channel annihilations into leptons (£ = p. r) 
through a Z' boson exchange and (b) f-channel annihilation into a pair of Z' bosons. 


The leading order DM annihilation cross sections are given by 

mj + 2 ml 


(™)w-+u - 


g ,4 Q% 




m 


27 t — 4 m |) 2 + 


1 - -4 + 0(0, 


m: 


{ av ) ipip->z'z 1 — 


mi-m^ 
4^ (m|, - 4 mfy ‘ 


m z' 


mz 


+ 0(v 2 ), 


where £ = p, r, p t , and z/ r . The decay width of the Z' boson is given by 

T z , ~ y 


9' 


,/2 


t=fi,T,Vp,V T 


12irmz' 


(m|, + 2m|) W1 


4 771 ^ 

e 6 (■ mz' — 2 me) 


m z> 


2 / 0/2 




4m3, 


127rm^ 

where 0 is the unit step function 


f M' + 2m l ) f / 1 - - 2m vO , 




(4) 

(5) 

( 6 ) 


B. Relic abundance 

Taking the DM relic density 0.11 < Ddm^ 2 < 0.13 [32], we found the preferred parameter 
space for -0 dark matter in mz> — g' plane for Q'^ = 2 in Figure |2j The plots for other values 
of Q[ p are also given later. The ballpark range is 1 < mz'[GeV] < 500 and 0.001 < g' < l.o 
as a reasonable choice within the perturbative regime. Notably, the dip structure appears 
around ~ 2= 20 GeV due to the resonance in the s-channel annihilation into leptons 
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FIG. 2. Preferred parameter regions in the mz> — g' plane for Qt = 2. In the narrow red band, the 
relic density of DM ^ is in the range of 0.11 < O^m^ 2 < 0.13. In the blue band, the annihilation 
cross section into /i + p~ and t + t~ satisfies + r + r - ~ (0.95 — 1.49) x 10^ 26 cm 3 /s, which 

is required to fit the GeV gamma-ray excess. 

mediated by the Z' gauge boson. In calculating the thermal average of DM annihilation cross 
section for the relic abundance, we take the the non-negligible effect of DM kinetic energy 
near the resonance pole, mz> — 2= 20 GeV as explained in Ref. |40j . The result is shown 
in Fig. [2] 


C. Fermi-LAT GeV excess 

We conduct the fit of our model, /i + /i~ : r + r” = 1:1, to the observed spectrum of the 
GC GeV y-ray excess. Our best fit is obtained for (crv)^^ + ~ 1.22 x 10~ 26 cm 3 /s 

with \ 2 — 19.22. Our result can be compared with the results in Ref. [Hj, where the best 
fits are (av) ~ 0.86 x 10~ 26 cm 3 /s with y 2 = 10.21 for the democratic leptonic final state 
and (av) ~ 1.42 x 10 _26 cm 3 /s with y 2 = 14.22 for branching ratios of : t + t~ = 

2 : 1. Accepting the y 2 < 29.6, which provides a p-value larger than 10 -3 for 10 degrees 
of freedom (i.e. 11 data points and one fitting parameter (av)), we found the preferred 
interval (av) = (0.95 — 1.49) x 10 -26 cm 3 /s. In Fig. [3j we plot our best fit as a red line 
with an interval corresponding to a p-value of 10 -3 . The data points are presented by 
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m#= 10 GeV, -» g + g & r + r 



E y [GeV] 


FIG. 3. Fits to the GC GeV 7 -ray excess for 10 GeV DM annihilating into g + and r + r~ with 
branching ratios of /i + /r~ : t + t~ = 1:1. The best fit is obtained with (av) ~ 1.22 X 10 _ 26 cm 3 /s, 
which is plotted as a red line. Upper and lower fits corresponding to a p-value greater than 10 -3 
are presented as purple dashed and dotted curves, respectively. The black points with blue error 
bars are the data points extracted from Ref. |14| . 


black dots and their error bars are represented by blue lines. As can be seen in Fig. [2j 
the plot for the right relic abundance reproduces successful GeV excess in the GC as was 
originally observed in [SI © and also in [Etj for leptonic annihilations taking into account 
the contributions by the ICS and bremsstrahlung with the annihilation cross section of 
(av) ~ (0.95 — 1.49) x 10~ 26 cm 3 /s for the preferred mass range near 10 GeV. In Fig. [2J for 
the case with Q), = 2, the parameter space fitting both the relic abundance and the GC GeV 
excess lies in three regions (m^'[GeV], g') = (< 9.6, 0.027), (19.7 — 20.3, 0.006 — 0.0012), 
and (30 — 42, 0.028 — 0.056). The parameter space would be slightly changed with different 
values of Q^. e.g. for = 0.1, (19.7 - 20.3, 0.026 - 0.0054) and (27 - 44, 0.10 - 0.27). 
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III. EXPERIMENTAL CONSTRAINTS FOR THE PREFERRED PARAMETER 


SPACE 

We now check whether the preferred parameter space mz> ~ 0(10 —100) GeV and g' < 1 
is still available after taking the relevant experimental constraints from the processes poten¬ 
tially induced by the gauged lepton number interactions: (g — 2) M , r decay, neutrino trident 
production, loop-induced DM-nucleon scattering and leptonically interacting Z' searches at 
colliders. 


A. (g - 2) m 


The gauged lepton number interaction leads to corrections to the muon anomalous mag¬ 
netic moment = (g — 2) M through a Feynman diagram shown in Fig. |d](a). The one-loop 
contribution to (g — 2) M is given by [4T1 [43] 


a 9' 2 

/\n - - L_ 

M 127r 2 rn%, ’ 


where we assume that mz> which is indeed valid with 

perimentally measured value and the SM prediction of (g — 2) 


rriz' ~ 10 GeV. 
M are respectively 


(7) 

The ex¬ 
given as 


(11659209.1 ±6.3) x 10" 10 , 

(8) 

(11659180.3 ±4.9) x 10“ 10 . 

(9) 


Thus, there exists discrepancy between the experimental value and the SM prediction: 

Aa /t = a® xp - a® M = (28.8 ± 8.0) x 10' 10 . (10) 


Even though the difference may be a sign of new physics but, more conservatively, we would 
set an upper bound on the size of the new contribution given in Eq. ([T]) and find the 2 a 
bound line in Fig. |8j 


B. t decay 

The gauged lepton number interaction may be seen in the leptonic decay of tau through 
the box diagrams such as Fig.Jd^b), which could make the branching fraction, Br(r —» /iu r z/ /t ), 






FIG. 4. Feynman diagrams that give a correction to (a) (g — 2)^ and (b) r —>• ^ decay. 


larger than what the SM predicts. It is interesting to notice that the measured value of the 
tau decay branching fraction to /iz/ r z7 /t is indeed slightly larger than what the SM predicted: 
the measured values for the branching ratio, Br(r —y yUZ/ r I7 M ), and the life time of tau from 
the PDG [33] are 


Br(r ->• ^ t ^)|pdg = (17.41 ± 0.04)%, 

r T | PDG = (290.3 ± 0.5) x 10~ 15 s , 


( 11 ) 

( 12 ) 


which has more than 2cr deviation from the SM prediction 


Br(r ->• nu T u^) 
Br(r ->• /in T z%) SM 


~ 1 + A with A = (7.0 ± 3.0) x 10“ 3 . 


(13) 


From the box diagrams with the Z' mediation, the deviation A could be evaluated 


A = 


3 g' 2 log(m^/?n| / ) 
47 t 2 1 — rn 2 z , / Tiil v 


(14) 


Interestingly, the sign of A from the U(1 ) L ^_ Lr interaction is consistent with that required 
by the difference between the experimental value and the SM prediction, Eq. (13). In Fig. [8j 
the upper region of an orange curve is excluded by the r —* nv T V^ decay limit at the 20 - 
level. 
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l/j 
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FIG. 5. (a) The leading order diagram with a Z l exchange contributing to neutrino trident pro¬ 
duction. (b) Diagram of dominant direct detection process. 


C. Neutrino trident production 


Neutrino trident production, , has been observed by several neutrino 

beam experiments such as CHARM-II [37] and CCFR 


OCHARM-II 

°SM 

CCCFR 


CSM 


= 1.58 ±0.57, 


= 0.82 ±0.28 


(15) 

(16) 


The measured cross sections are consistent with the SM prediction so that stringently con¬ 
strain our model. In the SM, the neutrino trident production is induced by a n + n~ pair 
production from the scattering of a muon-neutrino in the Coulomb field of a target nucleus 
[46] \W'\ . In our model, the leading order correction is coming from the contribution of Z’ 
boson shown in Fig. [5}[a) that interferes with the SM contribution from similar diagrams 
with a W/Z boson exchange instead of the Z’. In our analysis, we use the exclusion limit 
(95% C.L.) obtained from the CCFR data in Ref. [39] which is shown as a light cyan-shaded 
region with the cyan dot-dashed curves in Fig. [8] 


D. Dark matter direct detection 

DM direct detection experiments search for the recoil energy of nucleus by DM scattering 
off nucleus. In this model, DM does not directly couple to quarks at tree-level. However, 
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one-loop suppressed scattering processes such as the one shown in Fig. |5](b) can still provide 
a sizable DM-nucleus scattering cross section in spite of the loop suppression factor. The 
one-loop suppressed DM-nucleus scattering cross section is given by 


O’fN ~ 


1-1% r (oiemZ\ (m\ 

g l< 


-I 2 


(17) 


where A = mz> / (g' y Q'%) is the effective cut-off scale, hn = m^m^/ (mjv + m%) is the 
DM-nucleus reduced mass, and Z is the atomic number, i.e. the EM charge of the target 
nucleus. Note that Eq. © originally has a log dependence on the renormalization scale 
due to the fermion loop. However, such log dependences from //- and r-loops cancel each 
other out thanks to the relative sign difference between //- and r-loop induced diagrams. 
In order to directly compare the DM-nucleus cross section with experimental bounds, we 
convert Eq. (|l7|) into the DM-nucleon cross section using the following relation: 


_ J_/4_ 

Oim — A2 ,,2 (T ' I P N 
Ahv 


(18) 


where A is the atomic mass number of the target nucleus and p n is the DM-nucleon reduced 
mass. For m^ ~ 10 GeV, the most stringent direct detection bound is currently provided 
by the LUX experiment |51j . The LUX limit is shown as a purple dashed line in Fig. [8] 


E. Searches for Z —> \t at the LHC and LEP 

The LHC results also provide constraints on the gauged lepton number interactions 
through the lepton productions. A single Z' production is allowed at tree-level at hadron 
colliders such as the LHC in pp —>• i + i~ Z' where the Z' boson is radiated from a lepton in 
the Drell-Yan process as shown in Fig. [6] even though Z l interaction is lepton-specific. The 
produced Z' boson subsequently decays either to a pair of charged-leptons, neutrinos or DM 
particles: 

Z'^p + p~, T + T~, Vfp h V’V’ i (19) 

if kinematically allowed. These processes can be probed by detecting either one charged- 
lepton pair plus missing Et events or two charged-lepton pairs, i.e. 4£, at the LHC. In 
this work, we focus on the 4L signals due to its clean and distinctive signature. If rriz' 
(m T , mp), the branching ratios of the Z' become 

Br(Z' el) = Br (Z 1 ->• v t v e ) = Br (Z' -)• ^)!Q% . (20) 
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t ,u e /ip 


FIG. 6. Feynman diagram for a Z' boson production process at a hadron collider. The Z l boson 
is radiated from a lepton, and then decays into a pair of leptons or DM’s. 


The leptophilic Z' can be detected at the LHC in four charged-lepton final states. The 
dominant SM backgrounds for this process are 


pp trz trtr , 

( 21 ) 

pp — >■ zz — >■ trtr. 

( 22 ) 


Four charged-lepton (4£) production at the Z resonance has been already measured by 
ATLAS [52] and CMS |53| collaborations at the LHC. Three final states have been well 
observed: pp —> 4e, 2e2/i, 4p. We consider only the four muon final state since in our 
scenario the Z' does not couple to electrons. In this analysis, we use the following selection 
cuts which is used in the ATLAS analysis [52J: 

• Pt^ > 4 GeV and \rj\ < 2.7 for individual muons, 

• Separation of muons: AA /J/t > 0.1, 

• Invariant mass of a muon pair: M m + m - > 5 GeV, 

• Invariant mass of four muons: 80 GeV < < 100 GeV. 

In Fig. [7j we present the Z' production cross sections through the pp —>■ p + Z' process for 

g' — 0.1 at the 8 and 14 TeV LHC which is obtained using MadGraph 

This Z' can be produced at tree-level at lepton colliders such as LEP through the similar 
process as shown in Fig. [6] just by replacing qq with e + e~ since the gauged 11(1)^-^ boson 
also has no direct coupling to e + and e~. The potential constraint from LEP for rriz' < mz 
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FIG. 7. Z' production cross section at the 8 TeV and 14 TeV LHC, through pp —> g g + Z'. We 
have set g' = 0.1. 

has been well studied in Ref. [42]J through the process, e + e“ —> g + g~ Z'. Despite much 
smaller total integrated luminosity, the limit from LEP is comparable to that from the 8 
TeV LHC due to much cleaner signals. In Fig. [8j we present the LEP limit on Z' from 
Ref. [42] as a dark-gray shaded region with the black long-dashed curves. 


F. Summary of experimental constraints on the U(1 )l m -L t model 

In Fig. |8j we collectively depict all the relevant constraints to the gauged lepton number 
interaction in rriz' — g' plane. 

• The limit from (g — 2)^: We plot the 2a limit from (g — 2)^ as a green solid line 
in the mz' — g' plane for representative choices Q= 2, 1, 0.5, and 0.1. The upper 
region of the green line is constrained by the current measurements of (g — 2)^. 

• The limit from r —>• decay: The upper region of the orange curve is excluded 

by the r —>■ /i/y r z7 /t decay limit at the 2a level. 

• The limit from Neutrino trident production: The 95% C.L. exclusion limit is 
shown as a light cyan-shaded region with a cyan dot-dashed curve. 
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FIG. 8 . Allowed parameter space of the U(1 )l^~l t charged dark matter model in the mz> — g' 
plane for four representative DM charges Qf = 2, 1, 0.5, and 0.1 (from top-left to bottom-right), 
respectively. We present the regions satisfying the DM relic abundance 0.11 < Ddm h 2 < 0.13 
and the annihilation cross section + T + T - ~ (0.95 — 1.49) x 10~ 26 cm 3 /s required to fit 

the GC GeV excess as red and blue bands. The upper regions of green, orange, cyan dot-dashed, 
purple dashed, gray dotted, and black long-dashed curves are constrained by (g — 2) M , r decay, 
neutrino trident production, LUX, LHC, and LEP, respectively. In this analysis, DM mass is fixed 
as m^p = 10 GeV. 

• The limit from dark matter direct detection: The LUX limit is plotted as purple 
dashed lines for four representative values Q[f, = 2, 1, 0.5, and 0.1. 

• LHC Z — Ai limit: The light gray-shaded region with the gray dotted curve is 
excluded by measurements of the Z —> 4/i at the LHC [19], [50]. The Z —» 4/i searches 
at the LHC strongly constrain the parameter space of rriz' ~ 5 — 40 GeV since the Ai 
production has been measured at the Z resonance and the selection cuts of Pt,^ > 4 
GeV and M^+^~ > 5 GeV are used. 
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• LEP Z —> At limit: Dark grad-shaded region with the black long-dashed curve is 
excluded by measurements of the Z —> A/j, at LEP [32] . 

For Q[ p < 1, considerable parameter space has already been ruled out by neutrino trident 
production and Z —y 4/x observations at the 8 TeV LHC and LEP, except the region around 
the resonance of mz> ~ 2mq. In near future, for larger Qt > 1 most of preferred parameter 
space will be verified by DM direct detection experiments such as XENON1T. The region 
around the resonance will be complementarily proved by Z —> 4/x searches at the 14 TeV 
LHC. 

IV. CONCLUSION 

In this work, we have explored a leptophilic DM model with the gauged U(1 )lu-l t symme¬ 
try in the light of the Fermi-LAT GeV gamma-ray excess. With this simple leptophilic DM 
model, we can simultaneously explain the observed DM relic abundance and the Fermi-LAT 
GeV excess. Our leptophilic Z' DM model additionally contributes to the muon (g — 2), tau 
decay process, and neutrino trident production. In particular, neutrino trident production 
measurements provide the most stringent constraint to the DM model in most of the param¬ 
eter space. Despite the absence of direct couplings with quarks, this model can be strongly 
constrained by DM direct detection bounds through the loop-suppressed process. For DM 
with a large charge under the U(1 Q[i, > 2, the current LUX direct search limit is 

comparable or stronger than the neutrino trident production limit. The U(1 )l i 1 -l t gauge 
boson can be produced through the radiation process from Drell-Yan leptons, which has been 
constrained by Z —y 4/i searches at the LHC and LEP, especially for rriz> ~ 5 — 40 GeV. 
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